Infrared Studies of the Onset of Conductivity in Ultra-Thin Pb Films 
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In this paper we report the first experimental measurement of the infrared conductivity of ultra- 
thin quenched-condensed Pb films. For dc sheet resistances such that cut <C 1 the ac conductance 
increases with frequency but is in disagreement with the predictions of weak localization. We 
attribute this behavior to the effects of an inhomogeneous granular structure of these films, which is 
manifested at the very small probing scale of infrared measurements. Our data are consistent with 
predictions of two-dimensional percolation theory. 
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Transport measurements in ultra-thin films have been 
a subject of active interest over many years Q. These 
systems, consisting of a thin layer of metal deposited onto 
a substrate held at LHe temperatures, provide a rela- 
tively simple way to study the interplay between local- 
ization, electron-electron interactions, and superconduc- 
tivity in disordered quasi-2D metals. These experiments 
are in quantitative agreement with predictions of local- 
ization theory Q combined with the effects of diffusion- 
enhanced electron-electron interactions S. The reason 
why these theories, developed for homogeneous materi- 
als, work so well in the case of granular, inhomogeneous 
films is that the length scale at which electrons lose phase 
coherence in these measurements is usually much larger 
than the characteristic size of inhomogeneities (grains, 
percolation clusters, etc.) of the film. 

Yet another way to modify the length scale is by sim- 
ply changing the probing frequency. This gives rise to 
frequency dependence of the ac conductivity in the re- 
gion ujt <C 1, where the Drude theory predicts a plateau. 
While the experimental data on the frequency depen- 
dence of conductivity are virtually non-existent for ultra- 
thin quenched-condensed films, they abound for thicker, 
more granular films, deposited onto a warm substrate [Q. 
In ac conductivity the frequency itself defines a charac- 
teristic dephasing length scale = ^D/uj, where D is 
a diffusion coefficient. In the frequency range where L u 
is smaller than other dephasing length scales, it enters 
into all localization and interaction formulas and gives 
rise to frequency-dependent quantum corrections to the 
conductivity. However, these quantum effects constitute 
only one source of frequency dependence of the conduc- 
tivity. In the region where the material is strongly inho- 
mogeneous on the scale of , the frequency dependence 
of conductivity is dominated by purely classical effects 
due to charge dynamics on a network of capacitively and 
resistively-coupled clusters of grains. The effective way 
to describe these effects theoretically is provided by the 
framework of percolation theory H. In this theoreti- 



cal approach the ac conductivity is shown to increase 
with frequency. Indeed, since capacitive coupling be- 
tween grains is proportional to frequency, grains become 
more and more connected as the frequency is increased. 
It is also known that purely quantum effects such as lo- 
calization and diffusion-enhanced interaction corrections 
become profoundly modified on length scales where the 
material can no longer be treated as homogeneous [||. 
The presence of inhomogeneities is known to change the 
functional dependence of these corrections on the phase 
coherence length, while leaving the absolute magnitude 
of the effect virtually unchanged. 

In this letter we report the first measurement of con- 
ductivity at infrared frequencies in ultra-thin films. Films 
used in this experiment were made in situ by evaporating 
Pb onto a Si(lll) (sets 1 and 2) and glass (set 3) sub- 
strates, mounted in an optical cryostat, held at 10 K. Ag 
tabs, pre-deposited onto the substrate, were used to mon- 
itor the dc resistance of the film. Infrared transmission 
measurements from 500 to 5000 cm -1 (set 1), and 2000 to 
8000 cm -1 (sets 2 and 3) were made using a Bruker 113v 
spectrometer at the new high-brightness U12IR beamline 
at the BNL's National Synchrotron Light Source. The 
substrates were covered with a 5 A thick layer of Ge to 
promote two-dimensional thin-film growth, rather than 
the agglomeration of the deposited Pb in larger grains. 
Even when using this traditional method, which is known 
to promote the growth of more homogeneous films, where 
the continuity has been measured at near a monolayer of 
deposited metal, we only saw the beginnings of conti- 
nuity near 8 — 20 A of metal with consistently thinner 
films using the glass substrates. Films were evaporated 
at pressures ranging from the low 10 -8 to the mid 10 -9 
Torr range. The transmission spectra were obtained af- 
ter successive in-situ Pb depositions. The dc resistances 
in set 1 on Si range from 64MO/D at 17.4 A average 
thickness to 543 0,/D at 70 A. The 70 A sample was then 
annealed twice, first to 80 K, and then to 300 K. As a 
result its resistance at 10 K became 166 O/n after the 
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first annealing, and 100 f2/D after the second annealing. 
Films from set 2 (also on Si) are similar to set 1: we have 
observed R n = 20Mfi/D at 18 A and R n = 1000 fi/D 
at 88 A. Finally, films from set 3, deposited on a Ge- 
coated glass substrate, range from 13 to 200 A, while R a 
changes between 5.6 MSI and 22.8 SI. 

The transmission coefficient of a film deposited on the 
substrate, measured relative to the transmission of the 
substrate itself, is related to real and imaginary parts of 
the sheet conductance of the film as m 



V ' [l + Z ^H/(n+l)]2 + (Z 0( 7SM/(n+l)) 2 ' 

(1) 

Here Z$ = 377 Vt is the impedance of free space, n is the 
index of refraction of the substrate, equal to n$i = 3.315 
for silicon and uq = 1.44 for glass, and a' n (uo) (some- 
times called G), ct"(uj) are the real and imaginary parts 
of the sheet conductance of the film. Almost every- 
where in our experiments a' D (u>), cr'aioS) <C (n + 1)/Zq. 
In this case the contribution of the imaginary part 
of conductance to the transmission coefficient is neg- 
ligible and Eq. (@) can be approximately replaced by 
T(w) ~ [1 + Z a' D (oj)/(n + l)] -2 . Even for our thickest 
films, where cr'^(cu) w (n + \)/Z$, the error in calculating 
a' a (to) in this way is less than 10% over our frequency 
range. Throughout the manuscript we will use this ap- 
proximation to extract the real part of the sheet conduc- 
tance of the film from its transmission coefficient. Only 
for our thickest films will we use Eq. ([!]) to derive param- 
eters of Drude fits. 

In Fig. 1 we plot the frequency-dependent conduc- 
tance, extracted from the transmission data for the films 
from the set 3 with the help of the above approximation 
to Eq. (|l|) . The seven thickest films from this set exhibit 
a characteristic Drude falloff at high frequencies. For the 
rest of the films the conductivity systematically increases 
with frequency throughout our frequency range. The in- 
set in Fig. 1 shows the average ac conductance as well as 
the dc sheet conductance for set 3 as a function of thick- 
ness. Note the curves start to significantly deviate from 
each other at around 50 A. 

In order to fit the conductance of our thickest films 
with the Drude formula, one needs to use the untrun- 
catcd Eq. (|l|) for the transmission coefficient. Insert- 
ing the Drude expression for the sheet conductance 
an(u>) = &d/{^ — ioJT) directly into Eq. [l| one gets 
T(w)/[l-T(o;)] = (l+w 2 r 2 )/[( ( j D / ( 7o) 2 +2 ( j D / ( To],where 
(To = {n+l)/Zo. Therefore, the transmission data, which 
are consistent with the Drude formula can be fitted with 
a straight line, when T(u>)/[1 — T(w)] is plotted as a func- 
tion of lu 2 . In Fig. 2 we plot the 7 thickest films from the 
set 3 in this way. The knowledge of the average thick- 
ness of our films along with parameters of the Drude 
formula enables us to calculate the plasma frequencies in 



our films. They are shown in the inset of Fig. 2 as a func- 
tion of 1/<7d — the dc sheet resistance in Drude formula, 
which itself was extracted from our Drude fit. These re- 
sults are in excellent agreement with the experimentally 
determined lead plasma frequency of lu p — 59 400 cm -1 
B . In the remainder of the manuscript we discuss pos- 
sible interpretations of the increase of the conductance 
with frequency, which we observe in our thinner films. 

One mechanism which is known to cause a frequency 
dependence of conductivity within a Drude plateau 
(lot <C 1) is purely quantum mechanical in origin. The 
conductivity is known to be reduced due to increased 
back scattering of phase-coherent electrons (so called 
weak localization (WL) as well as diffusion enhanced 
electron-electron interactions (EEI) ||. The magnitude 
of this reduction depends on the length scale over which 
an electron maintains its phase coherence. In the ab- 
sence of external magnetic field or ac electric field this 
length is determined by the temperature. It is given by 
the inelastic scattering length li n (T) = y Dn n (T) for 
WL, and the thermal coherence length Lt = y/hD/kT 
for EEI. Here D is the diffusion coefficient of the elec- 
tron related to its dc conductivity by the Einstein for- 
mula a = e 2 (dN/dfj,)E F D, and n n (T) is the tempera- 
ture dependent inelastic scattering (dephasing) time. In 
the presence of the ac electric field the diffusive motion 
of an electron is restricted to a spatial region of size 
L u = y/D/u). If this length scale turns out to be shorter 
than the corresponding dc length scale, it is L u which 
enters in all WL and EEI formulas. 

The question of effective dimensionality of the quasi- 
2D sample is decided by comparing L u to the film thick- 
ness d. The frequency dependent WL corrections to the 
sheet conductance of the film are given by A<j 2D (u) — 
^j-lriLJT in the 2D limit (d < LJ) ||, and Aerg D (w) = 

j^j-dy/% in the 3D limit (d > L u ) ||. At the lower 
end of our frequency range u> = 500 cm" 1 for a realistic 
value of D = 5 cm 2 /s we can estimate L u < 20 A < d. 
Therefore, for our films one should use the formulas of 
three-dimensional localization theory. The frequency- 
dependent sheet conductance in most of our films is con- 
sistent with the y/uj dependence of 3D WL. However, we 
believe that in order to explain the frequency dependence 
of our experimental data one needs to look for yet another 
mechanism, supplementing that due to weak localization 
and electron-electron interactions. The problems with as- 
cribing the observed frequency dependence of conductiv- 
ity solely to WL and EEI effects are: (i) the dependence 
of the slope of the conductivity vs y/U on the thickness of 
the film and the dc sheet conductance, which determines 
the diffusion coefficient D, does not agree with predic- 
tions of the 3D localization, (ii) the weak localization 
theory is only supposed to work in the limit where its 
corrections are much smaller than the dc conductivity. 
In our experimental data we don't see any change of be- 
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havior as the corrections to conductivity become bigger 
than the dc conductivity. In fact the \fu> fit works very 
well and gives roughly the same slope even for films with 
dc sheet resistance of ss 100 kS7, while the ac sheet resis- 
tance is only as 1 kf2. Furthermore, the 3D-localization 
theory predicts that the t/uj dependence of weak local- 
ization theory should be replaced by u)( d ~ 1 '' d — cj 1 / 3 
dependence at or near the 3D metal-insulator transition 
p0| . In our experimental data we see no evidence for 
such a crossover. 

There exists yet another, purely classical effect that 
gives rise to the frequency dependence of the conduc- 
tivity. It is relevant in strongly inhomogeneous, gran- 
ular films. There is ample experimental evidence that 
even ultra-thin quenched-condensed films have a micro- 
scopic granular structure 11 1^]. In order to describe the 
ac response of a film with such a granular microstruc- 
ture one needs to know the geometry and conductivity 
of individual grains as well as the resistive and capac- 
itive couplings between grains. The disorder, which is 
inevitably present in the placement of individual grains, 
makes this problem even more complicated. However, 
there exist two very successful approaches to the analyti- 
cal treatment of such systems. One of them, known as the 
effective-medium theory [ fL3|] , can be viewed as a mean- 
field version of a more refined approach, based on scal- 
ing near the percolation transition. The EMT takes into 
account only concentrations of metallic grains and the 
voids between the grains, disregarding any spatial corre- 
lations. A more refined approach takes into account the 
geometrical properties of the mixture of metallic grains 
and voids. The insulator-to-metal transition in this ap- 
proach is nothing else but the percolation transition, in 
which metallic grains first form a macroscopic connected 
path at a certain critical average thickness d c of the film. 
The dc conductivity above the transition point scales as 
{d-dcf, where t = 1.3 in 2D and t = 1.9 in 3D @. Just 
below the percolation transition the dielectric constant of 
the medium diverges as e(d) ~ [d c — d)~ s , where s = 1.3 
in 2D and s = 0.7 in 3D. The diverging dielectric constant 
is manifested as the imaginary part of the ac conductivity 
o~(lu) ~ — iuj(f c — f)~ s . In general the complex ac conduc- 
tivity of the metal-dielectric (void) mixture close to the 
percolation transition is known J^] to have the following 
scaling form: 



a(uj, d) = \d— d c \ t F±(—iuj\d — d, 



!-(*+«) 



(2) 



Here F+ (x) and F _ (x) are scaling functions above and 
below the transition point correspondingly. Note that 
this scaling form correctly reproduces the scaling of 
the dc conductivity above the transition and the di- 
vergence of the dielectric constant below the transition 
provided that F + (x) = f} 0) + F^x + f} 2) x 2 + . . ., 
while F-(x) — F^x + F^x 2 + . . . One should men- 
tion that the predictions of the EMT can also be writ- 



ten in this scaling form with mean-field values of the 
exponents t = s = 1, and scaling functions F±(x) = 
(V-D 2 + 4(D - l)x±D)/(2(D- 1)), where D is the spa- 
tial dimension. 

Since the metallic grains in our films form not more 
than two layers, our data should be interpreted in terms 
of the two-dimensional percolation theory. In two dimen- 
sions t = s = 1.3 ||, and according to Eq. (|^) the ac con- 
ductivity precisely at the transition point d = d c is given 
by a(u,d c ) = A(iu>/oj ) t/( - t+s) = A(iu> /w ) 1/2 ■ This pre- 
diction is in agreement with our experimental data. In 
Fig. 3 we attempt the rescaling of our data according to 
Eq. (|J). The critical thickness d c is determined as the 
point where the ac conductivity divided by ^/Zj is fre- 
quency independent. Of course, the experimental uncer- 
tainty in our data points does not allow us to determine 
which exponents t and s provide the best data collapse. 
However, as we can see from Fig. 3 our data are consis- 
tent with the scaling form of the 2D percolation theory. 
Finally, we use Fig. 3 to estimate basic parameters such 
as typical resistance R of an individual grain and typ- 
ical capacitance C between nearest neighboring grains. 
From the limiting value of a(uj, d)(d c /\d— dd) 1 ' 3 at small 
values of the scaling variable x = u(d c /\d — d c |) 2 ' 6 for 
d > d c , one estimates the conductance of an individual 
grain to be of order of R ~ 1000 £1. In the simplest RC 
model, where the fraction of the bonds of the square lat- 
tice are occupied by resistors of resistance R, while the 
rest of the bonds are capacitors with capacitance C, the 
ac conductivity exactly at the percolation threshold is 
given by A/R^luRC) 1 ^ 2 , where A is a constant of or- 
der of one. Therefore, the slope da/d^/uj in our system 
should be of the same order of magnitude as V 'RC ' / R. 
This gives C ~ 2.6 x 10~ 19 F, which is in agreement with 
a very rough estimate of the capacitance between two is- 
lands 200 A x 200 A x 30 A separated by a vacuum gap of 
some 20 A, giving C~2.7xl0 19 F. This order of mag- 
nitude estimate confirms the importance of taking into 
account inter-island capacitive coupling when one inter- 
prets the ac conductivity measured in our experiment. 
Indeed, R = 100017, and C = 3 x 10~ 19 F define a char- 
acteristic frequency 1/RC ~ 17000 cm^ 1 comparable to 
our frequency range. 

In summary, we have measured the conductivity 
of ultra-thin Pb films in the frequency range 500 to 
8000 cm -1 . The evolution of a(uj) with DC sheet re- 
sistance is consistent with classical two-dimensional per- 
colation theory in this range. At lower probing frequen- 
cies, where becomes larger than the scale of inhomo- 
geneities in these films, we expect that the effects of weak 
localization will become more prevalent. 

We have benefited from fruitful discussions with 
P.B. Allen, A.M. Goldman, V.J. Emery, V.N. Muthuku- 
mar, Y. Imry, Z. Ovadyahu, and M. Pollak. The work at 
Brookhaven was supported by the U.S. Department of 
Energy, Division of Materials Sciences, under contract 



3 



no. DE-AC02-98CH10886. Support from NSF grants 
DMR-9875980 (D.N.B.) and DMR-9725037 (B.N.) is ac- 
knowledged. Research undertaken at NSLS is supported 
by the U.S. DOE, Divisions of Materials and Chemical 
Sciences. 



FIG. 3. The rescaling of conductivity data in set 3 accord- 
ing to predictions of the percolation theory. The value for d c 
is 34 A. 
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FIG. 1. Sheet conductance vs frequency for set 3. The 
dashed lines plotted between 3000 and 4000 cm -1 (where the 
glass substrate is opaque) are a guide to the eye. The inset 
shows the inverse average conductance in our frequency range 
(filled symbols) and the dc sheet resistance (open symbols) as 
a function of the film thickness. 



FIG. 2. T/(l-T) plotted vs. co 2 for the 7 thickest films 
from set 3; The solid lines are Drude Model fits, as described 
in the text. The inset shows the plasma frequency for these 
films (open symbols) and three quenched films from sets 1 
and 3 (filled symbols) as a function of the dc sheet resistance, 
with the flat line representing the plasma frequency of bulk 
lead obtained from Ref. 8. 
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